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The maximum power point tracking (MPPT) has been a popular terminology 
among those researchers who deal with one of the most available renewable 
energy resources called solar power. Many researchers have evolved and 
proposed a lot of techniques for tracking maximum power. But each 
technique has its pros and cons. Some techniques suffer from complexity as 
far as implementation is concerned, whereas other techniques lack accuracy. 
In current research work, it has been investigated the MPPT using Newton 
Raphson (NR) method, which has not yet received a good attention by 
researchers. After observing its limitation, this method is modified (i.e., 
abbreviated as modified Newton Raphson method (MNRM)) to make it 
suitable for extracting the maximum power from PV module. The feasibility 
and precision of this method depends upon accurate measurement of 
temperature and irradiation. By using MNRM, the presumption of initial 
value is close to the voltage corresponding to maximum power, so it 
leads faster converging to solution through a few rounds of 
iterations. In order to validate it, a MATLAB/Simulink model is developed 
and simulated. The proposed method is incorporated in PV module-fed buck 
converter so as to explore superior performance. 
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1. INTRODUCTION 


A photovoltaics PV cell is smallest unit of solar power generating unit and its basic characteristic 
(I~V) is mostly non-linear in nature. The connection of PV cells in series and parallel is known as PV module 
and the grouping PV modules in series and parallel is known as PV array. The modeling of PV cell is 
presented by Salmi et al. [1]. Since DC-DC converter is major component in solar power unit, so the 
comprehensive review of various DC-DC converters is detailed and analyzed [2]. The maximum power point 
tracking MPPT is common phenomenon among researchers in the field of solar technologies and the various 
techniques of MPPT are focused and compared [3]. In this literature, the comparison includes the positive 
aspects of each technique and its limitation. Also, it has been discussed about their performances and 
complexities about their implementation. Anand and Fernandes [4] have elaborated on implementation and 
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power management of solar photovoltaic system. In literature [5]—[7], the researchers have presented the 
design, simulation and implementation of grid-tied solar power controller integrated with micro-inverter for 
power supply. The design and simulation for grid connected solar PV system connected with battery backup 
is presented [8]. For improvement of power quality [9], it is implemented the solar PV system associated 
with battery storage system based micro-grid operation. The researchers [10], [11] have successfully 
implemented solar PV-based multi-level inverter for power conditioning in distribution system and induction 
motor drives for water pumping system respectively. It is also presented analog, multiplier MPPT controller 
and optimal controller for solar PV based dc power supply for improvement of performances [12], [13]. The 
comparative study of maximum power point tracking under variable size perturbations is presented [14] and 
findings are placed. The concept of fuzzy logic controller [15] is applied in case of photovoltaic system for 
micro-grid system. A novel algorithm of particle swarm optimization (PSO) is extended solar PV system for 
better performance [16]. Nasir et al. [17] have critically investigated for rural electrification in developing 
nations, The application of single phase to three phase inverter and smart inverter are extended solar PV 
system [18], [19]. Here it has been focused the method of connectivity of solar inverters. Ogudo and 
Umenne [20] has presented design for buck-boost converter interfaced with solar PV based power supply. 
The output voltage of this converter is non-inverting type. A novel adaptive maximum power controller [21] 
is proposed for photovoltaic system and the parameters of controllers are automatically tuned to meet the 
performance requirement of PV system. On the other hand, it is proposed a model predictive controller [22]. 
This is exclusively for bidirectional DC/DC and AC/DC converters. This indicates that both variable output 
DC/AC and regenerative power feedback/power reversal are possible. Above all, the whole literature review 
comprises various controllers and modified Newton Raphson method MPPT methods [23]—[30]. In current 
research work, the MPPT technique is developed using MNRM that leads to faster convergence to the 
solution, which is later treated as tracking value. This MNRM is incorporated in PV module-fed buck 
converter and the whole system is simulated with MATLAB/ Simulink. 


2. ANALYSIS OF PV MODULE 

The solar PV array is formed by connecting the PV modules in series and parallel combinations to 
acquire the required current and voltage rating. The sum of the individual power rating of the module is equal 
to the rating of the PV array [31].The basic goal of full power point tracking is to interpret the voltage and 
current from the solar panel, measure the power, and then display the power to its limit. Multi-stage DC-DC 
converter manages MPPT, and PWM controls grid current from Inverter. The inverter circuit provides all the 
activities in a single phase that involve MPPT and grid current power [32]—[37]. 

The mathematical model of PV system is described by (1) the photo current. 


lon = [Ise + K,(T — 298)] G/1000 (1) 
The saturation current is (2). 


q Ego(1/Tn - 1/T) 


I, = Iys (T/T) exp [ nK ] (2) 

The reverse saturation current (3). 
Isc 
mS a 
e\N Ns KT) _4 
The output current I = Inn — la — Ish 
(V+I Rs) (V+IRs) 

or I = lpn — Nolo expC Ea oa (4) 

The power is P = VI 
(V+I Rs) (V+IRs) 
PSLE Vin ee a a e (5) 


Where 

Ipn= photo-current (A) 

Is= Short-circuit current (A) 

ki = Short circuit coefficient of the cell(at nominal temperature 
25 degrees centigrade and at G= 1000 W/m’) = 0.0032 
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T : Operating temperature (degree kelvin) Tn : Nominal temperature (degree kelvin) = 298 
G: : Solar Irradiation (W/m?) q_: Electron charge (C) = 1.6e-19 coulomb 
Voc: Open-circuit voltage (V) n : Ideality factor of the diode =1.3 

K : Boltzmann’s constant (J/K) = 1.38 e-23 Ego: Band gap energy of the semiconductor (eV) =1.1 

Ns : Number of cells in connected in series Np : Number of PV modules connected in parallel 

Rs : Series resistance (Q) Ryn: Shunt resistance (Q) 

The parameters of a specific module are given as follows and considered in forthcoming simulation work. 
Isc =8.21A, ki = 0.0032, Voc =32.09, N5=54. Np=1, R= 0.221 Q, Ren = 415.405 Q, G <=1000 W/m’. 

In Figure 1, the characteristic of power, current and load conductance (Y) of PV unit with respect to 
its terminal voltage (Vpy) without MPPT are shown. The load conductance (Y=I,,/Vpy) goes on decreasing 
from maximum towards zero with respect to gradual increase in terminal voltage of the circuit from zero to 
open-circuit condition(V<).The load resistance is reciprocal of load conductance. So higher load conductance 
leads to lower load resistance and vice-versa. With reference to peak power, there corresponds to current 
(Imp), voltage (Vip) and conductance (Ymp). SO Ymp is the load conductance that corresponds to maximum 
power. Between peak power to open-circuit voltage, the power falls drastically, as the conductance tends 
towards zero. Similarly, from short-circuit position (i.e., v=0) to peak power, there is gradual change in 
increase of power. The I~V characteristic is complex in nature and upon close attention of graph of current 
and power with respect to voltage, it is noticed that the current fall is significantly small between v=0 to V mp 
and it is approximated as current source. Beyond Vmp, there is drastically fall in current, but the voltage 
difference between v=Vmmp and Voc is small, so characteristic beyond v= Vmp is approximately treated as 
constant voltage characteristic. 

The MPPT concept has been popular among researchers since more than a decade to explore the 
possibility of maximizing power at different load. The researchers have evolved various techniques for 
tracking the point of maximum power and the review paper [3] pertaining to this concept is illustrated. All 
these techniques lead to determination of tracking maximum power point so as to maximize power for 
particular load or to utilize the power as per researcher’s choice. In current research work, the tracking of 
maximum power point is accomplished using MNRM and hence to accomplish desired task, a control 
strategy is developed which is explained in next section. 

The buck converter is known as a prominent member among family of DC-DC converter having 
action similar to DC-DC chopper. This is also utilized now-a-days as a tool in extracting power from solar 
power unit. With controlled switching of the buck converter, output power can be controlled for desired 
value. It is assumed that input and output voltage remain constant over the cycle. 


ag 


Figure 1. Curves of current, power and conductance with respect to voltage of PV unit at irradiation 
(G=1000 W/m?) and temperature (T=30 °C), scale: Ypv-(1:10) , Ipv (1:10), Ppv (1:1) 


3. CONTROL STRATEGY 
3.1. Determination of Vmp by modified Newton Raphson method 

In the determination of Vmp (i.e., voltage corresponding to maximum power) is necessary to know 
the boundary between constant current and voltage region. Though various techniques have been evolved by 
researchers earlier, but in this case, the MNRM technique is utilized to solve this problem. 


Xx) 
" (Xx) 


(6) 


Xk+1 = Xk- 
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Where Xx and Xx+ is the kn and (k+1)tn iteration value of the f(X) = 0; It is needed the substitution of initial 
value of the solution Xo in above equation. The iteration is continued till the difference |X 441 — X:{ is very 
small, which is order less than le-3. For the solution of non-linear equation f(X), there arise two major 
constraints. One is the substitution of proper initial value, without which there may occur divergence of the 
solution of f(X) instead of convergence. The second constraint is that many a times, during iteration of above 
equation, it may so happen that the derivative of f(Xx) in denominator of above equation may be zero. So, 
this may lead to non-arrival /failure to the solution. 

Now these two problems can be attended using following ways. For a fixed PV module, the open- 
circuit voltage, Voc is highest and fixed with its highest irradiation and lowest temp. The main objective 
considering lowest temp is that the Voc would increase with decrease in temperature. Secondly, the voltage 
Vinp Which corresponds to maximum power lies very close to open-circuit voltage Voc than to it short-circuit 
position (V=0). Also, the variation in temperature has a least significant impact than variation in irradiation. 
So based upon above facts, the Voc corresponding to highest irradiation and normal temperature can be 
determined from mathematical model of PV module or its experimental setup. The accuracy of determining 
Voc is dependent upon the accurate measurement of parameters for PV module. So, this Voc can be used as 
initial value in solving NR Method and as this value is very close to V mp and as a result, there will be faster 
convergence of solution for its non-linear equation. This is the way to attend the Ist problem. The solution 
for 2nd problem is that in order to overcome the problem for the derivative of f(Xx) getting zero in the 
denominator of (6), it is added with a tolerance (£), for which the equation in is modified as (7). 


Xen = Xk pa (7) 
Where ‘e’ a small value that varies between le-4 to le-6. 

The initial value of variable for solution is the Voc (i.e., the maximum possible open-circuit voltage 
Voc) and it is treated as Xo. The inclusion of attending problems with above two ways in NRM is referred as 
MNRM and it would act as full proof for solution. It has been frequently observed that if the substitution of 
initial value Xo is close to real solution, then it occurs faster converging to solution. 


3.2. Development of mathematical model non-linear equation form for solution of Vinp 

The determination of Vmp that corresponds to maximum power is determined by computing 
dP/dV=0 using (5). 

dP/dV = d(IV)/dV = I + V(dI/dV) = 0 (8) 


Substituting value of ‘I’ from (4) into (8) and it is simplified as (9). 


(1+cV)exp(cV)=k (9) 
wherek = pi tlo isn y raphi 
T No Nolo 
a Gia ER ED gd 
nKN, T nKN, T 


Since Ip, is function of T and G and these two variables are determined by proper sensor, then constants‘k’ and 
‘c’are determined and this leads to solution of non-linear equation. In (9) is simplified as (10). 


f(V) = (14+ cV) exp(cV) —-k=0 (10) 


In (10) is the non-linear equation that undergoes MNRM for solution of getting Vmp that corresponds to 
maximum power. 


3.3. Flow-chart for implementation of control strategy 

The flow-chart for above control strategy is shown in Figure 2. Initially, the measurement of 
temperature and irradiation is to be measured and these two parameters are necessary to determine the Vmp, 
voltage corresponding to peak power, from MNRM. This voltage is assumed as reference voltage of buck 
converter. The buck converter is to be operated under close-loop control during constant voltage region (i.e., 
between v= Vinp and Voc) of I~V characteristic of PV unit. The load voltage and the current through filter 
inductance are considered as two feedbacks of close-loop control scheme. The error from comparison of 
reference voltage with actual one is processed through PI controller and then limiter to generate reference 
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current. This reference current is then compared with those of filter inductance and then processed through 
hysteresis current controller to generate gate pulses of switching device. 

For implementation of control strategy, the flow-chart shown in Figure 2 is transformed to 
MATLAB/Simulink/State flow. The complex logical thought process associated with MNRM is put into 
State flow block of Simulink. The Simulink model for whole control strategy is shown in Figure 3. 


Measure Temperature (T) 
and irradiation (G) of 
PV unit 


Development of the non-linear 
equation from mathematical 
modelling of PV unit. 
Apply MNRM to determine Vmp 


Take Vmp as reference for Buck 
converter under close-loop control 
scheme with voltage feedback 


Measure the error between 
the voltage across load 
and reference 


Process this error through PI 

A controller to derive current reference. 

This reference is compared with actual current 
in filter inductor and error is passed through 

hysteresis controller to generate gate pulses 
for switch in Buck Converter 


Figure 2. Flow-chart of control strategy 
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Figure 3. Close-loop control circuit of PV-fed Buck converter usingmodified Newton Raphson method 


4. RESULTS AND DISCUSSION 

The simulation results from conventional PV module without any close-loop is presented in Table 1. 
It shows that the variation of current, power and load conductance with respect to variation in voltage. The 
results of simulation from close-loop control of PV-module fed buck converter is given in Table 2. It shows 
the load voltage and power with respect to variation in load conductance. 

The results of power from Tables | and 2 are plotted in Figure 4. It shows the power profile with 
respect to load conductance under open-loop and close-loop. Similarly, from the results of Table 2, it was 
plotted a graph for load voltage and power with variation in load conductance in Figure 5. From the 
observation, it shows that the load voltage remains constant during power from its zero to its maximum value 
and the load conductance (x-axis) corresponding to maximum power is Ymp. But after the peak power, the 
close-loop controller fails to oblige. The observation is critically investigated and findings are noted as follows. 


Table 1. Power and load voltage of PV module with respect to variation in load conductance 
(without MPPTand buck-converter) (i.e. at T=30 °C, G=1000 W/m?, Vref=26.87 V) 
S1. No Load conductance (Y)mho Power without MPPT and buck converter(W) Load voltage (V) 


1 0.010 9.58 30.97 
2 0.020 19.00 30.83 
3 0.040 37.30 30.54 
4 0.05 46.2 30.4 
5 0.066 60.6 30.15 
6 0.100 87.85 29.64 
7 0.200 155.0 27.84 
8 0.294 187.3 25.24 
9 0.33 185 23.59 
10 0.5 133.7 16.35 
11 1 67.24 8.2 


Table 2. Power and load voltage of PV modulewith respect to variation in load conductance 
(with MPPT and buck-converter) (i.e. at T=30 °C, G=1000 W/m?, Vref=26.87 V) 
Sl. No Load conductance (Y) mho Power with proposed method (W) Load voltage (V) 


1 0.010 7.18 26.8 
2 0.020 14.05 26.74 
3 0.040 28.09 26.53 
4 0.05 35.11 26.46 
5 0.066 45 26.1 
6 0.100 71.87 25.97 
7 0.200 121 24.6 
8 0.294 171 24.17 
9 0.33 174.0 22.85 
10 0.5 132.84 16.3 
11 1 66.91 8.18 
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Figure 4. Output power of PV module (Ppv) under Figure 5. Voltage and power profile with respect to 
open loop and (Pload) under close-loop (i.e. output variation in load conductance of PV module-fed 
voltage of PV module-fed buck converter) with buck converter under close-loop operation 
respect to variation in load conductance 


The close-loop control of buck converter with input having PV unit is stable under during V=V mp to 
Voc of I~V characteristic of PV module, which can be characterized as constant voltage (CV) 
characteristic/region and it is equivalent to Y=0 to Ymp of Figure 5. But in case of constant current 
characteristics (i.e., from V=0 to Vmp of I~V module), the constant current causes filter inductor to saturate 
early and thus fails to boost the voltage across load. This region is same as falling power with variation in 
load conductance in Figure 5. So, the former one, which is the constant voltage region of PV module, is 
considered to interface buck converter with PV unit for stability. Also, the load conductance to be considered 
in this region must lies in the range (i.e., between Y=0 and Y=Y mp) so as to facilitate the satisfactory 
operation of close-loop control. The design of filter inductance and capacitance is based upon switching 
frequency and load condition. As the switching frequency is small and variable due to its operation under 
hysteresis current controller, so upper range of switching frequency is considered to design filter parameters. 
The design values of filter inductance and capacitance across load are considered as 10 uH and 500 uF 
respectively. The waveforms of current through filter inductance, filter capacitor, load resistance and voltage 
across load are shown in Figure 6 for two cycles under steady-state condition during close-loop operation 
with load resistance of 10 Ohm. The current through the filter inductance (Iis) is the summation of currents 
through filter capacitor (Ics) and load resistance (Io). 


0 
0.038 0.0381 0.0381 0.0381 0.0382 
Time ---— >(sec) 


Figure 6. Waveforms of currents through filter inductance, filter capacitor, load resistance and voltage across 
load under steady-state condition during close-loop operation 


5. CONCLUSION 

A new concept of tracking maximum power using MNRM is successfully implemented. Here the 
issue of assumption of initial condition is sorted out using open circuit voltage of PV module, which is very 
close to the voltage (Vmp) corresponding to maximum power. As the assumption of initial value in MNRM is 
close to the voltage corresponding to maximum power, hence it leads faster converging to solution through a 
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few iterations. But the specialty of this technique is to monitor continuously the measurement of temperature 
and irradiation of PV unit and the performance of this technique is dependent upon accuracy of sensors 
tracking above two parameters (i.e., temperature and irradiation). Also, the other feature of this proposed 
method is that the close-loop operation is implemented for the load conductance varying from zero to Y mp or 
during constant voltage region of PV module that corresponds to maximum power. Beyond this range, the 
controller fails to take any action, which is the limitation of the proposed scheme. This scheme is applicable 
mostly for standalone dc load, where the matter of stabilizing accurate output voltage is not a matter of great 
concer. 
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